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Abstract 

We calculate the specific heat of the ideal gas obeying the generalized ex- 
clusion statistics (GES) in the continuum model and the tight binding model 
numerically. In the continuum model of 3-d space, the specific heat increases 
with statistical parameter at low temperature whereas it decreases with sta- 
tistical parameter at high temperature. We find that the critical temperature 
normalized by /ij (Fermi energy) is 0.290. The specific heat of 2-d space was 
known to be independent of g in the continuum model, but it varies with g 
drastically in the tight-binding model. From its unique behavior, identifica- 
tion of GES particles will be possible from the specific heat. 
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1 Introduction 

The many-particle wave function is symmetric under an exchange of two identical par- 
ticles for Boson and is antisymmetric for Fermion. The dimension of the Hilbert space 
is determined by counting the wave functions constrained by the symmetry. Intermediate 
statistics are under the intensive study recently. In one and two-dimensional space, arbitrary 
phases can appear in many-particle wave functions under an exchange of two identical parti- 
cles, which are called anyon . In contrast with the anyon, Haldane proposed generalized 
statistics ten years ago without specific reference to spatial dimensions [H. It is called the 
generalized exclusion statistics (GES). In his proposal, the dimension of the Hilbert space 
d, and the particle number A^, are connected by 

Ad = -gAN (1) 

where Ad is the change of the available single-particle states number and AA^ is the change 
of the particle number for identical particle system, g is a statistical parameter. We call 
particles obeying the GES as g-on. The state-counting of many particle states obeying the 
GES is proposed later by Wu, which corresponds to Boson with g = and Fermion with 
(7 = 1^. This counting rule seems to be exact in the limit of a large number of states [Q, 
and the distribution function of the ideal g-on gas has been derived. 

Because the statistics is the fundamental property of nature, it is important to know if 
there are particles obeying generalized statistics. For this purpose, physical properties of the 
ideal g-on gas should be fully understood. Thermodynamic properties of the ideal g-on gas 
are those properties and have been studied by several people. The mean occupation number 
was derived in the context of the one- dimensional solvable model ||^,|^ and the lowest Landau 
level anyon model |p. Wu derived the entropy and free energy and so on and Nayak et 
al. showed the duality The Sommerfeld expansion was applied PJ9|-[TT|. Other thermal 



properties were studied in [|T2|. In particular, the specific heat of the ideal g-on gas, which is 
the most fundamental observable quantity which may reflect the statistics, has been studied 
analytically [^^j-Q in a simple model where the single particle energy is proportional to p^. 



The g-ons which have more general energy dispersions have not been studied so far. In the 
present work, we find the specific heat of the ideal g-on gas in several spatial dimensions of 
wide energy band system and narrow energy band system based on numerical calculations. 

First, we study the wide energy band system where the single particle energy is pro- 
portional to p^. In three-dimensional space, we find that the specific heat increases with 
statistical parameter at low temperature region whereas it decreases with statistical parame- 
ter at high temperature region. The critical temperature normalized by (Fermi energy) is 
found to be 0.290, at which the specific heat is independent of g. In two-dimensional space, 
the specific heats are independent of the statistical parameter. The analytical proof of the 
(^-independence of the specific heat was given in [§,0. The origin of the ^f- independence is 



due to the constant density of state (DOS) for ideal gas in two-dimensional space. 

Next, we study the ideal g-on gas on lattice space model, which has the narrow energy 
band, in two-dimensional space. The DOS is not constant in this model. This case would be 
realized in the fractional quantum Hall state (FQHS). The quasiparticle appearing in FQHS 
IB|^T5[ is known to be an anyonic soliton with a finite spatial extension. Two quasiparticles 



cannot overlap each other and have a minimum distance. So it would be reasonable to 
associate these quasiparticles with the particles on a lattice and to describe with the tight- 
binding model. It was shown, in fact, that anyons in the lowest Landau level behave as ^f-ons 
with a varying DOS ||^. Thermodynamic properties of these gas have not been studied. 
Thus we study the specific heat of ideal g-on gas in the tight-binding model. By numerical 
calculations, we find that the specific heat of this model is totally different from the previous 
cases and varies with a statistical parameter of the GES. 

The present paper is organized in the following way. In section 2, we review the GES 
and give the analytical expression of the distribution functions. In section 3 numerical 
calculations of the specific heats of the ideal g-on gas are given in three and two-dimensional 
space. In section 4, the specific heat of the tight-binding model of g-on gas is calculated 
numerically. Summary and discussion are given in section 5. In Appendix A, we give some 
useful relations. In Appendix B, we give the exact solution for the g-on distribution function 



of g = where is a natural number. Improvement of the Sommerfeld expansion and the 
duahty of the coefficient are given in Appendix C. 

2 Distribution function of g-on 

For Boson and Fermion, the number of quantum states of identical particles occupying 
G states is given by 

- N\{G-1)\ ' 
'''' = NKG-NV. 

respectively. Wu proposed the interpolation formula implying the GES as 

» Nl[G - gN - {1 - g)]\ ^> 

with g = corresponding to Boson and g = 1 Fermion. This is the state-counting rule for 
g-on. 

Using Eq. (||), the distribution function reads 

^^'^ = y{e^(^-^))+g 
where the function y satisfies the functional equation 

yiCy{l + yiOY~' = C, C = e^(^-^). (6) 

From Eq. (|]) y{() = ( ior g = 1, and y{() = ( — 1 for g = 0, respectively. We can easily 
obtain other exact solutions of Eq. (j^) for g = 2,3,4, by elementary calculation. 

Moreover, the other solutions for g = where is an arbitrary natural number (A^ > 1) 
are given by 

ri = ( A \ \ 1 F^r of J- J- 2. ^-2 _J, 2_ N~2 -iV^e^(^-P) _ ll /'7^ 

y \^ J^q j [A'-l-' W-iViV' ATI JV'"' JV ' iV-l ' AT-i >■■■ JV-I ' (jv-l)JV-l / -^J \' ) 

where x = f , P = f and iFn{oii, ■ ■ ■ ,ai, (3i, ■ ■ ■ , z) is the hypergeometric function. The 
temperature t is normalized by /i/ = j). The derivation of y is given in Appendix B. 



The positivity of ( leads to the generahzed exclusion principle 

/(e) < (8) 

At zero temperature, we get C = and y{C) = for e — //^ < 0, and ( = oo and y{C) = oo 
for e — fig > 0, where fig is a pseudo- Fermi energy for arbitrary g cases. Then the value of 
the average occupation number at zero temperature is 

I i (e-/i„ < 0), 
[O (e-/i, >0). 

Eq. (P) means that one particle occupies ^ states at zero temperature. 

3 The specific heat of the ideal g-on gas in 3-d and 2-d space 

We numerically compute the specific heat of the ideal g-on gas in 3-d and 2-d space in 
this section. We assume that the particle of any GES has the same mass and the spectrum, 

2 

and neglect the spin degree of freedom. In (i-dimensional space, the DOS for the system 

d — - 

with particle number, A^, is D(e) = q(d)e2~^ where q{d) = ^fJ^f ^ and fif is Fermi energy 
given by Eq. (0) in Appendix A. The chemical potential fi{T) is determined by the particle 
number = D{e)f{f]{e — fi))de. Using the /i(T), we can calculate the average energy E 
and the specific heat C from next relations E = D{e) e /(/3(e — fi))de and C = —ksP"^^ 
where j3 = -j^. Converting the integral variables etoy for N and E gives the next relations 

d = '\L — yfa+i) 



and 



N^f 2 Ja y{y + l) 

a corresponds to the value of y ioi e = 0, which satisfies — = g\oga + {\ — g) log (a + 1). 
For arbitrary ^'-on except Boson, we can solve Eq. (p^ ) with respect to a at fixed temperature 
t numerically and obtain the average energy from Eq. (0). 
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d = 3 

For three-dimensional space, the average energy E{t) is shown in Fig. |I]. At low temper- 
ature, the average energy increases with g, which reflects the generalized exclusion principle. 
At high temperature, all curves have the same slope and they go to the classical limit. 

The specific heat versus temperature is shown in Fig. 0. The specific heat changes 
continuously and has one fixed point at t = 0.290. The entropy increases with g below this 
temperature, but decreases with g above this temperature. This refiects the generalized 
exclusion principle. For g = case (Boson), infinite particle numbers in the ground state 
cause the Bose-Einstein condensation. For g case, no condensation occurs at low 



temperature as mentioned in [|T0[ . The behavior of the specific heat is different for different 
statistics, so in principle we are able to identify the GES from the specific heat. 

Some arguments [^P-p!T[l, in which the specific heat of the ideal g-on gas is expressed 
as a power series of the temperature (Sommerfeld expansion for Fermion), suggest that the 
coefficient of the first power of the temperature increases with g in 3-d space (we improve 
the Sommerfeld expansion and show the duality of the coefficient in Appendix C). Our 
calculation is not based on the Sommerfeld expansion, and it is possible to check the validity 
of the Sommerfeld expansion. The deviation between the linear part in the temperature and 
numerical result is given in Fig. |^. At low temperature, the deviation is small, so the 
Sommerfeld expansion is good. However, the deviation decreases with g. For the Boson's 
case below the Boson's critical temperature it behaves as t^^^. So it may be natural to see 
that the specific heat for very small g differs from the linear temperature dependence. Hence 
the Sommerfeld expansion is not good for small g. 

d=2 

For two-dimensional space, the average energy E{t) is plotted in Fig. ^. At zero temper- 
ature, E reflects the generalized exclusion principle and depends on g. At high temperature 
it approaches to the classical limit. The speciflc heat is shown in Fig. ^. Figure ^ shows that 
the speciflc heat is independent of the statistical parameter g at arbitrary temperature. The 
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analytical proof of the (^-independence of the specific heat has been given in PJTT[] , using the 



e-independence of the DOS (-D(e) = ^ for 2-d space). 

4 Specific heat of the tight-binding model of the ideal g-on gas in 2-d space 

The quasiparticle of Laughlin's theory of the FQHS has a finite energy width and is ap- 
proximately described by the lattice model. Girvin et al. [|16| calculated the energy spectrum 



of the quasiparticle of the Laughlin's theory numerically, based on the single mode approx- 
imation analogous to the Feynman's theory of superfiuid helium. Their result implies that 
the energy spectrum of the excited state in FQHS seems to be explained by the tight-binding 
model. Laughlin's theory describes the FQHS {v means the filling factor of the FQHS) of 
v = \ beautifully. In this case, the quasiparticle has | charge |1T3|JT5[| , and satisfies fractional 



statistics of |. Hence there may be three quasiparticles in one available single-particle state; 
that is (7 = I statistics. It would be reasonable to express the quasiparticle in the Laughlin's 
theory by the tight-binding of ^^-on gas. We study the g-on in the tight-binding model, 
consequently. 

In the tight-binding model, the energy spectrum is e=-c{cos{k^b)+cos{kyb)) and the DOS of 
this model is given by 



where c is the hopping constant and b is the lattice spacing, and M(= ^) is the number of 
the lattice. K{x) is the complete elliptic integral of the first kind 

K{x) = rd<P (13) 
•^0 Jl-x'^sin (b 



The particle number and the average energy of this model is given by 



N= r deD{e)f{e) = — r dexiJl 

J-2c CTT^ J~2c \ V 



2cJ / yiO+g' 



E= deD{e)ef{e) = — de K\ J 1 ^ ^ 



-2c CTT^ 7- 2c \ V \2 cj j y{Q + g 
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It is very difficult to calculate the specific heat for arbitrary value of g. Hence we select 
special values of g, g = 1, |, |, We fix ^ = |, which corresponds to the half-ffiling state 
for Fermion and the chemical potential of Fermion is zero. We numerically calculate the 
temperature dependence of the average energy by using and obtain the specific heat. 

The result is given in Fig. ||. We cannot obtain the specific heat of g = ^,0 at low 
temperature range now owing to technical problem of numerical calculation. In the high 
temperature limit, the one-particle energy goes to zero. In this limit, the distribution of the 
low energy particle number is very small, because many particles have much of energy. The 
integral decreases at high temperature, since the integral region is restricted between —2c 
and 2c. As a result, the specific heat vanishes in the low temperature limit. On the other 
hand, the specific heat becomes maximum at one temperature value of order 1; that is, when 
the temperature goes to the order of the hopping constant (the temperature is normalized 
by the hopping coupling c), the fluctuation of energy is very large. From the DOS in Fig. 
1^, we read that the most of the contribution comes from the particle having the energy of 
the order of hopping constant. This explains why C{t) has a maximum qualitatively. C{t) 
becomes maximum at t = 0.45,0.56,0.61 for g = The maximum values are 0.49, 

0.71, 0.77 (Fig. y). The difference of the maximum values reflect each statistics. See Fig. |^. 

From our results, g = case, which will be realized in the ^ = \ FQHS, is different 
from Fermion case. Therefore, we will be able to observe the exotic statistics 5* = | by 
measuring the specific heat of the ^ = \ FQHS at the temperature for the order of the 
hopping constant. However, it would be difficult to observe the specific heat of 2-d electron 
system, since the thermal effect of the orthogonal directions to the plane affects the plane's 
thermal phenomena, in which the FQHS is realized. 

5 Summary and Discussion 

The continuum model in three-dimensional space shows that the specific heat depends on 
statistics. Therefore, we will be able to observe the signature of exotic statistics by measuring 
the specific heat. Moreover, we find that there is a critical temperature t = 0.290 by 
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numerical calculation. Above this temperature the specific heat increases with the statistical 
parameter but below this temperature the specific heat decreases. 

As was shown before, the continuum model in two-dimensional space shows that the 
specific heat does not depend on statistics. The (^-independence is caused by the constant 
DOS. In the tight-binding model, where the DOS is not constant, the specific heat of the 
ideal g-on gas depends on statistics even in two-dimensional space. Hence we are able to 
distinguish the GES by measuring the value of the specific heat especially from the peak of 
the specific heat. The temperature and the value of the specific heat in which the specific 
heat becomes maximum is obtained for several g values. It would be exciting to identify the 
exotic statistics. 

T.-H. A. is grateful to K. Ishikawa for suggesting these subjects and careful reading of this 
manuscript and also thank N. Maeda and J. Goryo for inspiring and critical discussions. It is a 
pleasure to thank H. Suzuki for very useful discussions. 

Appendix A: Some useful relations 
From Eq. (|^) y satisfies 

y{Cni + y{Oy-' = C (14) 

where ( = e^^'^~^\ Taking logarithm of it gives 

(3{e- ii)= g\ogy+{l- g)\og{y + l). (15) 
We convert the integral variable e into y. From Eq. (plSf ) we obtain 

* = ^?^*- 

The density of state D{e) of the ideal gas in ci-dimensional space is given by D{e) = 
^^r^y (^) where h is the Planck constant. The Fermi energy is determined by = 

J,^' deD{e) is 



where P = y {V = L'^). Especially, in three and two-dimensional space, the Fermi energy is 
given by 

Appendix B: Exact solution of the distribution function for g = jf 

We obtain the analytical solution of y in this Appendix In particular, in the case of 



g = where A is a natural number (except A = 1), y is represented by the hypergeometric 
function. 

In general, the i-th solution Zi of f{z) = 0, where z is a complex variables, is gained by 

where the contour of the integral is circle around Zi. From Eq. (^) in low temperature 
limit and e — /i < 0, we see that e~t'~ — > and then y ~ e~t^. So the solution of y in low 
temperature limit will behave as e~«^ and go to zero. Using this formula and analytical 
continuation of y, the solution is given by 

1 f , + + y{y+l)^-' 



y = — fdy — y (20) 

where the contour of the integral is circle around and x = |, p = j. In low temperature 
limit, Eq. (pOD is expanded by The residue gives the integral value. After short 

calculation, the solution y is given by 

Because the singularity of the solution depends on g, we consider the case of g = where 
A is natural number to avoid the singular expression of Eq. (0). By using the important 
relation of Gamma function T(z)T(i-z) = the solution is represented by a nonsingular 

form. The procedure of the calculation is as follows; 
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^ n=l r(-t+2) n\ 



E 



{ T{nN-n) 



oo / ^N{x-p)\n , iv-i _ 



r(n-^)r(n)r(n-i+A)...r(n-i+v: 



n=l 



n! 



(Ar_l)(jv-i)n 

TT^ 



(27r)^-v^]V^ 



r(n)r(n+^)r(n+^)---r( 



1 /AT-l g ^ ^ _ArA^e^(-"P) \ " r(n-i)r(n-i+ a) . . . r( 



1 I iV-1 



27rivV iv n! V (iv - ir- V r(n+^)r(n+^) ■ ■ ■ r(n+#5f 



1 



'TV- 1 



27rA^ V 



/(0,iV) 



1 / -A^^e^("-P) \" 7(n,iV) 



(iV>l) 



1 



The solution of y is 



N-2 . -N" e 



AT'JV'iV' JV 'JV-l'JV-1' JV-l' (jY-l)-'V-l 



1 



9 i i ^ 

V JV > JV ' JV ' 



JV-2 . 1 2 
JV 'JV-l'JV- 



JV-2 . 
■ JV-l ' 



_jv'V^JV(a:~p) ■ 



(JV-l) 



N-1 



(22) 



where > 1. In Fermion case (A^ = 1), we return to the first expression (^) and can get 
the Fermion distribution function. The hypergeometric function is represented by multiple 
integral, so by analytical continuation the present solution is valid in all temperature region. 
This solution perfectly corresponds to = | (A^ = 2) and 9 = 1 = 3) which is 
obtained exactly in Eq. (|^). We numerically calculate the low temperature behavior of the 
3-d specific heat by using this solution. 



Appendix C: Duality for the coefficient of the Sommerfeld expansion 

We improve the Sommerfeld expansion and show the duality relation of the coefficient 



for g-on in this Appendix . The Sommerfeld expansion is the expansion in terms of small 
deviations from the step function 



X 



s-l 



dx = 

y + g JO 



1 \ x^~^ x^~^ 

-] dx + / dx + / ( 

q / Jp y + g Jo q 



-dx. 



We set 
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p = (1 -^)log2 

to express the upper and lower value of the integral of y explicitly. Then the expansion is 
as follows; 

ax = X \ ax + / ax + / ax 



y + g Jo \y + 9 9/ Jp V + 9 Jo 9 

= - fx-^^^dx + r ^dx + ^J- 

Jo 9[y + 9) Jp y + 9 9 s 

= / ip-uY'^—r-^ — -du+ f (p + vY~^^—dx + -^. 
Jp 9{y + 9) Jo y + g g s 

In the limit of t — {p ^ 00), we can approximate this integral as follows 

/ dx^- ip-uY-' , , du+ (p + vY-^ dx + -- 

Jo y + g Jo g{y + 9) Jo y + g 9 s 

^,T{s-n)T{n+i)^ ' g Jo y + 9 
n=o r(s-n)r(n+i) Jo y + g gs 

. _ £ (ti):B„ to) + c„ to)) + '-^ 

^0 r(s-n)r(n+i) V y s 

where —x + p = u, x— p = v in third line, and Bn{g) and C„(5f) are given by converting u, 
V into y 



B^)= ru-^du= f -J^{-g\ogy-{l-g)\og{y + l) + {l-g)\og2)\ 

Jo y + g Ji y + l 

/■oo 1 /■oo (Iqi 

C„(.)^/ dv^ y^^^glogy+{l-g)\og{y + l)-{l-g)\og2f. 

Jo y + g Ji y{y + l) 
In Bn{g) and Cn(s), we convert y into 2; = and obtain 

5nW= r -^(-^logz + log(l-2) + (l-<7)log2r, 

Jl/2 1 — Z 

C„(.)= /' -(^log^-log(l-z)-(l-5)log2)". 

Jl/2 Z 

Moreover, in C„(fl) and Bn{g), we convert z into 2; = 1 — w and into ^ respectively, then 

dw 

'1/2 



Cn(s) = - / (^log(l-«;)-log«;-(l-^)log2r, 

Jl/2 1 — W 

^-W = - / 7— (-log^ + 5log(l-2)-(l-5)log2r. 

0" >/l/2 1 — 2; 
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So we get the duality relation 

Finally, the Sommerfeld expansion is given by 

Jo y + g ^ r(s-n)r(n+i) V ^ " ) g s 

The Fermion case (g^ = 1), the odd number of n in the sum does not cancel. 
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FIGURES 




FIG. 1. Temperature dependence of energy in 3-d space. The vertical axis is e = and 
the horizontal axis is temperature normalized by The broken line represents the Boson case 
(5 = 0). 
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FIG. 2. Temperature dependence of specific heat in 3-d space. The vertical axis is c = and 
the horizontal axis is temperature normalized hy fif. The dashed line represents the Boson case, 
and the dashed and dotted line represents the Fermion case. The critical temperature is 0.290. 
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FIG. 3. Comparison of the Sommerfeld expansion with the numerical calculation in 3-d space. 
The first term in the expansion is written by a solid line and numerical results are given by a 
dotted line. The deviation is large for small g. The first term in the expansion is = 
Note that the scales of t are different for each figures. 
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FIG. 4. Temperature dependence of energy in 2-d space. All lines are parallel. They differ by 
constant. From top to bottom, g = 1,^, |, 0.1, 0.01. 
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FIG. 5. Temperature dependence of the specific heat in 2-d space {g = 1, 3,0.1,0.01). Ah 
hnes overlap. 
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FIG. 6. The distribution function at t = 0.5. The vertical axis is the distribution function /(e) 
and the horizontal axis is one-particle energy normalized by the hopping constant (e = ^). 
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FIG. 7. The product of /(e) and d{e) at i = 0.5. The vertical axis is the product /(e) and 
e) = and the horizontal axis is one-particle energy normalized by the hopping constant. 
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FIG. 8. The specific heat of the tight-binding model of the ideal g-on gas. The vertical axis is 
E and the horizontal axis is temperature normalized by the hopping constant. 
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FIG. 9. The temperature where the specific heat becomes maximum and its value. 
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